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Abstract 

Due to quantum correlations in the C-odd and C-even D°D° pairs produced in the reactions 
e + e~ — ► D°D°(mr ) and e + e~ — > D D ^(ni: ), respectively, the time-integrated D°D° decay rates 
are sensitive to interference between amplitudes for indistinguishable final states. The size of this 
interference is governed by the relevant amplitude ratios and can include contributions from D°-D° 
mixing. We present a method for simultaneously measuring the magnitudes and phases of these 
amplitude ratios and searching for D°-D° mixing. We make use of fully- and partially-reconstructed 
pairs in both C eigenstates, and we estimate experimental sensitivities based on a plausible 
charm factory dataset. Similar analyses can be applied to coherent K°K°, B°B°, or B®B® pairs. 

PACS numbers: 14.40.Lb, 13.20.Fc, 12.15.Mm 
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I. INTRODUCTION 



Studies of the evolution of a K° or B° into the respective anti-particle, a K° or B° [l|, 
have guided the form and content of the Standard Model and permitted useful estimates of 
the masses of the charm @, Q and top quark 0, Q prior to their direct observation. Neutral 
flavor oscillation in the D meson system is highly suppressed within the Standard Model 
and, thus, with current experimental sensitivity, searches for D°-D° mixing constitute a 
search for new physics. In addition, improving constraints on charm mixing is important for 
elucidating the origin of CP violation in the bottom sector. 

The time evolution of the D°-D° system is described by the Schrodinger equation 

™ iJ\ fD°(t)\ 
M - -T 1 




2 J \D°(t) J ' 

where the M and T matrices are Hermitian, and CPT invariance requires Mu = M 22 = M 
and I\i = r 22 = T. The off-diagonal elements of these matrices describe the dispersive or 
long-distance and the absorptive or short-distance contributions to D°-D° mixing. We define 
the neutral D meson mass eigenstates to be 

|Di) = p\D°) + q\D ) (2) 
\D 2 ) = p\D°)-q\D ), (3) 

where \p\ 2 + \q\ 2 = 1, and, following Ref. 0, D 1 is the CP-odd state, and D 2 is the CP-even 
state, so that CP|P°) = —\D°). The corresponding eigenvalues of the Hamiltonian are 

Xi,2=M 1)2 ~r lt2 =(M~rj±^ [m 12 - 1 -v 12 ^ , (4) 

where M 12 and r 12 are the masses and decay widths, respectively, and 



Q 



M* — IF* 

M 12 2 1 12 (5) 



p \m X2 - |r 



12 



A D° can evolve into a D° through on-shell intermediate states, such as K + K~ with mass 
Mk+k- =M d o, or through off-shell intermediate states, such as those that might be present 
due to new physics. This evolution through the former (latter) states is parametrized by 
the dimensionless variables — iy (x). We adopt the conventional definitions of these mixing 
parameters: 

M 2 - M l 

% = f (6) 

r 2 — Ti 

y = (7) 

The mixing probability, Rm, is approximately (x 2 + y 2 )/2 0]. For hadronic flavored final 
states, the above time evolution is also governed by the relative magnitudes and phases be- 
tween Cabibbo-favored (CF) and doubly-Cabibbo-suppressed (DCS) amplitudes, generically 
denoted by r and S, respectively. 
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Standard-Model-based predictions for x and y, as well as a variety of non-Standard- 
Model expectations, span several orders of magnitude [8J. Several non-Standard Models 
predict |x| > 0.01. Contributions to x at this level could result from the presence of new 
particles with masses as high as 100-1000 TeV 0, E3]. The Standard Model short-distance 
contribution to x is determined by the box diagram in which two virtual quarks and two 
virtual W bosons are exchanged. Next-to-leading order calculations show that the short 
distance contributions to x and y are expected to be comparable [ll[ . Long-distance effects 
are expected to be larger but are difficult to estimate. It is likely that x and y contribute 
similarly to mixing in the Standard Model. 

Measurement of the phase 7/^3 of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing 
matrix [l2| is chal leng ing and may eventually be limited by experimental constraints on 
charm mixing [l3|, 14]. Several methods have been proposed using B T — > DK T decays: 
the Gronau-London-Wyler (GLW) [TEj method, where the D decays to CP eigenstates; the 



Atwood-Dunietz-Soni (ADS) |16| method, where the D decays to flavor eigenstates; and the 
Dalitz-plot method [171 ] . where the D decays to a three-body final state. Uncertainties due 
to D decays contribute to each of these methods. The CLEO-c physics program [l8[ includes 
a variety of measurements that will improve the determination of 7/^3 from the 5-factory 



experiments, BABAR and Belle |19l . [20j. The pertinent components of this program are: 
improved constraints on charm mixing amplitudes (important for GLW), first measurement 
of the relative strong phase Sxn between D° and D° decay to K + tt~ (important for ADS), 
and studies of charm Dalitz plots tagged by hadronic flavor or CP eigenstates. The total 
number of charm mesons accumulated at CLEO-c will be much smaller than the samples 
already accumulated by the 5-factories. However, the quantum correlation of the D°D° 
system near threshold provides a unique laboratory in which to study charm. 

The parameters x and y can be measured in a variety of ways. The most precise constraints 
are obtained by exploiting the time-dependence of D decays [7]. Previous attempts to 
measure x and ^include: measurements of the wrong-sign semileptonic branching ratio 



D° -> K + £-j?i 
D° -> K+K~ 
sign D° — > K~ 
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24 



and 7r + 7r" 



which is sensitive to Rm', decay rates to CP eigenstates 
which are sensitive to y; the wrong- 
which is sensitive to x' 2 = 



26, 27 



29J30 



31 



33. 34 



~7r hadronic branching ratio 
(y sin 5k-k + x cos <5rv) 2 and y' = y rcos 8k-k — x sin Sktt] wrong-sign D° 
D° 
6 



K° S 7T' 



7T 



[37J, 



K+n-ir [35] and 
which determines 



K + 7t~7t + tt~ [36j decays; and the decay rate of D° - 
- from a Dalitz-plot analysis and measures x and y. 
Time-dependent analyses are not feasible at CLEO-c; however, previous authors have 
found that the quantum-coherent D°D° state (38, 39, 40] provides time-integrated sensitivity, 



K °sl 



through the interference between amplitudes for indistinguishable final states, to y at (9(1%) 
and to cos^xtt at 0(0.1) in 
extend the work of Refs. 



1 fb" 1 of data at the ^(3770) @, [fl Sj. In this paper 



we 

and develop a method for simultaneously measuring 
x, y, r, and 5. Unlike the proposed measurements of Ref. [6|, we do not rely on external 
estimates of the relevant D° branching fractions. Our method is a modified version of 
the double tagging technique originally developed to measure D branching fractions at the 
V>(3770) 0, EI, S3 . We make use of rates for exclusive D°D° combinations, where both D 
final states are specified (known as double tags or DT), as well as inclusive rates, where either 
the D° or D° is identified and the other neutral D decays generically (known as single tags 
or ST). Although we estimate that CLEO-c will not have sufficient sensitivity to observe 
Standard Model charm mixing (see Section II Vj) . it should be able to achieve a precision 
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comparable to current experimental results. The analysis presented in this paper can also 
be applied to coherent K°K°, B°B°, and B®B® systems, although with some additional 
complications. 



II. FORMALISM 

As in Refs. 0, Hl|, 42, 43], we consider the following categories of D° and D° final states: 



/ or /: hadronic states that can be reached from either D° or D° decay but that are 
not CP eigenstates. An example is K~tt + , which is produced via CF D° transitions or 
DCS D° transitions. We include in this category Cabibbo-suppressed (CS) transitions 
as well as self-conjugate final states of mixed CP, such as non-resonant K^tt 



7T 



• £ + or I : semileptonic or purely leptonic final states, which, in the absence of mixing, 
tag unambiguously the flavor of the parent D. 

• S + or S-: CP-even and CP-odd eigenstates, respectively. 

All D° decay modes can be treated uniformly if we enumerate charge-conjugate final 
states separately, indexed by j and J. For instance, K~tt + and K~£ + U£ are labeled by j and 
K + 7i~ and K + £~i?e by J. CP eigenstates appear in both lists because J = j. We define the 
mode-dependent amplitude ratio (j\D°) / (j\D°) = rje~ t5j , with 



T 3 



U\D°) 



(j\D°) 



( (J\D°) \ 



-5j = arg yjj^jj = Strong + <5wcak + 7T, (9) 

where the phase of ir results from our CP convention. Since 5 wea k i n the charm sector is trivial 
(0 or tt), 5j corresponds to either — Strong (if j is CF) or vr — Strong (if j is CS). Furthermore, if 
CP is conserved, then (j\D°) / (j\D°) = (j\D°) / (j\D°) . To resolve the ambiguity of whether 
to identify any given final state as j or J, we choose < Tj < 1. 

Since x and y have both been constrained to be less than 0(1%) [II], we generally keep 
terms only to leading order in x and y. We denote decay amplitudes by Aj = (j\D°) and 
Aj = (j\D°). In our phase convention, CP conjugate amplitudes are given by 

(10) 

(11) 
(12) 

We use a normalization in which A 2 - is the D° — > j branching fraction in the absence of 
mixing. If mixing is present, then the branching fractions for an isolated neutral D meson 





(f\D°) = - 


-(f\D°) 


A, = 


(£ + \D°) = 


-(r\D°) 


As ± = 


(S±\D°) = 


T(S ± \D°) 



produced in a D° or D° flavor eigenstate become 43 



Bj = B(D° —>■ j) m A 2 -{\ + rjfjj) (13) 
B, = B(D° - J) » A){r) + rrfj) = B 3 R 3 , (14) 
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TABLE I: Values of the amplitude ratio magnitudes rj and phases 5j, as well as uncorrelated 
branching fractions, for each final state category, to first order in x and y. 



j 


r 3 


Sj 


f 


r f 


Sf 


I 


r f 


Sf 


i+ 







s+ 


1 


TT 




1 






Aj(l + r f y f ) 
A 2 f {rj+r f y' f ) 
A 2 

V) 



A 2 l(l + y) 



where yj = y cos Sj + xsin5j, y', = y cos 5j — x sin 5j, and Rj = T(D° — > j)/T(D c 



The total D decay rate is unaffected by mixing, so 

E (4 + 4) = E 4 (l + r|) = E Bj (1 + Rj) = 1. 



(15) 



If the deviation of q/p from unity is parametrized by two small CP- violating parameters 
(magnitude and phase), then these parameters only appear in products with x and y; they 
can only modulate the strength of the mixing signal. Therefore, below, we assume q/p = 1 
and also that CP is conserved in the decay amplitudes (i.e., \(j\D°)\ = \(j\D°)\), which 
allows y to be expressed in terms of A 2 : 



E2^r,cos5 / = E4 + -E^ 



E 2A / r / cos5 /> 
/ 



(16) 



where we have accounted for the fact that S± modes are simultaneously labeled by j and j. 
Table [T] lists the values of rj and Sj for each final state category. 

As shown in Ref. 47J], a D°D° pair produced through a virtual photon in the reaction 
e+e" -> D°D° + rwy + nvr is in a C = (-l) m+1 state. Thus, at the ^(3770), where 
no additional fragmentation particles are produced, there is only C-odd, while at higher 
energies above D*D threshold, we can access both C eigenstates. The DT rates for final 
states j and k are given by jrj 41. 42. 43l 



r c -(j,k) = Q M \A(-\ 3 ,k) 



R 



M 



B ( -\j,k) 



T c+ (j,k) = Q'u AW(j,k) +R' M BW(j,k) + C™(j,k), 



(17) 



where 



A^(j,k) 
B (±) (j,k) 

CM(j,k) 



(j\D )(k\D )±(j\D°)(k\D ) 

P(j\D )(k\D )±V(j\D°)(k\D ) 
q p 

2^\A^*(j,k)B(+\j,k) ^ u 



IX 



(1-y 2 ) 2 (1 + x 2 ) 2 



(18) 
(19) 

(20) 
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Q't 



M 



R't, 



M 



i-y 2 i 
i + y 2 



X' 



[i-y 2 ) 2 
i + y 2 



+ 



X 



'1 + x 



2\2 



X 



[i-y 2 ) 2 

Using Equations M and M, we find 

l^ (±) 0',^)| 2 ^ |^? (±) (j, A;)| 2 
I^ (±) (J,A;)| 2 



(1 



„2\2 



A 2 A 2 



Qm -x 2 +y 2 
3R M - 



A 2 A 2 



\B {± \j,k)\ 2 
C^(j,k) « A)A\c% 
C^(j,k) 



1 + r 2 r\ ± r jrk v jk 
r) + r\ ± r jrk v± 



A 2 A 2 p- 



(23) 
(24) 



(25) 
(26) 
(27) 
(28) 



where Zj = 2 cos 5^, lo, = 2sm<5,,-, = (z.,-Zfc ± WjW k )/2, and 



2^2 



(1 + r 2 )r k z k + (1 + rDrjZj 



± 



[l + x 



2\2 



(1 - r 2 )r k w k + (1 - rl)rjWj 



(29) 

In Tables |TT| and IIII[ we give the D°D° branching fractions to DT final states for C- 
odd and C-even initial states, evaluated using the above formulae. If both D° and D° 
decay to the same final state, we divide the l^^l 2 and C*- 4 "** terms by 2. In Table [TTT| the 
entries with vanishing rate would be non-zero only the presence of both mixing and CP 
violation [42|, |43[ . If the D°D° decay were uncorrelated, these DT branching fractions would 
be B(j, k) = B(j, k) = BjB k (l + RjR k ) and B(j, k) = B{j, k) = BjB k (Rj + R k ), with a factor 
of 1/2 if j = k in the latter expression. 

The D°D° inclusive branching fraction to the ST final state j is obtained by summing all 
DT branching fractions containing j and is found to be the same for C-odd and C-even (to 
first order in x and y) and simply related to the isolated D° branching fractions: 



B(j,X)=-£[B(j,k)+B(j,k) 



A) 



1 + r • + rjZjy 



Ba + B T . 



(30) 



Table HVl shows these ST branching fractions evaluated for the three categories of final states. 
The total D°D° rate is obtained by summing either DT or ST rates: 

rw = E [r0'> k ) + *)] + E r (j, ^) = \ll \P(j, x ) + r (J> *)] ■ (3i) 

j,k>j j,k j 

Like the total rate for an isolated D°, the total D°D° rate is unaffected by mixing and 
quantum correlations. 

For the C-odd configuration, with only one mode of type /, the ST and DT rates depend 
on only four independent parameters: rf, Zf, y, and Rm- For the C-even configuration, there 



TABLE II: D D DT branching fractions for modes containing / or /, to leading order in x and y. 



f 
f 
f 

r 

£+ 
tr 

S+ 

S- 



C = -1 



A}R M 
A A 



l + rj(2 



l + ri(2 



4) 



AjAj, (r 2 f + r 



A)R M 



l + ri(2 



r f r f' v ff 



A 2 A 2 



f"f' i 1 + r / r /' " r f r f' v ff> 

A 2 A 2 r 2 
A 2 A 2 

A 2 f A 2 s+ [l + r f (r f + z f )] 
A 2 A 2 _ [l+r f {r f -z f )\ 



AjAj, (l + rjrj, 
AjAj, ( 



rj + rj, 



-z 2 ) + rj 

- TfTfiVjf, 

r f r f ,v+ f , 



A 2 A 2 
A 2 A 2 r 2 

J i f 

A 2 f A 2 s+ [l + r f (r f + z f )] 
A)A 2 _ [l + rf (r f -z f )} 



f 
f 

r 
r 

r 
s + 



C = +l 



4 4 



2A 4 

„2 



f r f 



r f + y' f + rjy f 



1 - rj(2 - zj) + rj + 4r f (y f + rjy' } )\ 2Ajr f (r f + y' f + r}y f 

A ) A )> ( r } + 4 + r f r f'*ff> + 2c 7/') A ) A )> i 1 + r ) r )' + r f r f' v ff> + 2c 
AjAj, (l + rjrj, + r f rf>vj f , + 2cj f , 



AjA 2 (rj + 2r f y>) 
jA 2 (l + 2r f y f 
*/)] 



AjAj, ( 



rj + rj, + TfTfiV 

\2 A2 



ff 



//' + 2C //' 



AjAi(l + 2r f y f ) 



A 2 A 2 

f t 



l-2y) 



[1 + 77(17 
[l+r / (r / + z / )](l + 2y) 



(r| + 2r /y ' f ) 



[l + r/(r/-^)](l-2y) 
[l + r / (r / + z / )](l + 2y) 



is one additional parameter, w/x, which appears via and y/. So, there is, in principle, 
sensitivity to x from the C-even configuration, although no information can be gained if Sf 
is or 7r. In addition, estimates of R~m and y can be combined to obtain x 2 . However, in 
all of these cases, it is impossible, without knowing the sign of 5f, to determine the sign 
of x. The mixing and amplitude ratio parameters can be isolated by forming ratios of DT 
rates and double ratios of ST rates to DT rates. Table |V] lists a selection of these ratios and 
functions thereof, evaluated to leading order in rj, x, and y. 



III. EFFECT ON BRANCHING FRACTION MEASUREMENTS 

If quantum correlations are ignored when using coherent D°D° pairs to measure D° 
branching fractions, then biases may result. For instance, if a measured branching fraction, 
denoted by B, is obtained by dividing reconstructed ST yields by the total number of D°D° 
pairs (A/ - ), then B differs from the desired branching fraction, B, by the factors given in 
Table ED 



Using a double tag technique pioneered by MARK III [44J, |45|, CLEO-c has recently mea- 
sured B(D° -»■ K-7T+), B(D° -> K-n + n ), and B(D°^> R-ti+ti+tt-) in a self-normalizing 
way (i.e., without knowledge of the luminosity or D°D° production cross section) [46-], us- 
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TABLE III: D°D° DT branching fractions for semileptonic modes and CP eigenstates, to leading 
order in x and y. 



t 



S- 



1+ 

£- 
S- 



A\R M 

4 4 
A t 

A 2 A 2 
A 2 A 2 



C = -1 



A^Rm 
A 2 A 2 







S + 
S- 
S' + 
S'_ 



3AjR M 
A\ 

A 2 A 2 + (l-2y) 
A 2 A 2 (1 + 22/) 
AfA 2 , + (l-2y) 
A 2 A 2 , (1 + 2y) 



C = +l 



3AjR M 

2A* s+ (l-2y) 
1 + 2y) 
A 2 A 2 l+ (l-2y) 4A 2 + A 2 /+ (l-2y) 
A 2 A 2 (1 + 2y) 



A 2 A 2 + (l-2y) 

A 2 A 2 
A e A s 



2^|_(l + 2y) 


4A 2 _A 2 , (l + 2y) 



TABLE IV: D D inclusive ST branching fractions, to leading order in x and y. 



j 


c = 


+1 and C = - 


1 


f 


A 2 


1 + r 2 + rfZfy 




I 




A 2 




s± 




A 2 s± (lTy) 





ing C-odd D°D° pairs from the ?K3770). Measured ST and DT yields and efficiencies are 
combined in a least-squares fit [49] to extract the branching fractions and TV. Quantum 
correlations were not explicitly accounted for in this analysis, but their effects were included 
in the systematic uncertainties. Only flavored final states were considered, and DCS contri- 
butions to the ST yields were removed, so the observed branching fractions are: 



r°-(/,f) 



B f (l 



r fyf- r f>yr - r f r r v ff> 



l-2r f y f + rj(2-zj) for /' = /. 



(32) 
(33) 



Similarly, the relationship between the observed M c , which is used to obtain the D°D° 
cross section, and the desired Af c ~ is 



Af' 



c _ T c -(f,X)T c -(f,X 

r c -(fJ') 



1 +r f y f + r f ,y r +r f r r v ff , 



(34) 



1 + 2^-^(2-4)1 for /' = /. (35) 
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TABLE V: Selected ratios of DT rates and double ratios of ST rates to DT rates, evaluated to 
leading order in r? , x, and y. Rates are represented by the notation = T(j, k) and Tj = T(j, X). 





C-odd 


C-even 


(1/4) • (T es+ T s jT es Ss + ~ rtS-Ts + /Tts + Ts-) 


y 


-y 


(T fi -/AT f )-(TsjTes.-r s+ /T es+ ) 


y 


-y 


(r ff -/Ar f ) • (r s JT fs _ - r s+ /r fs+ ) 


y + r f z f 


-(y + rfZf) 


(r/r s+5 _/4) • (i/r /5 _r 5+ - i/r /5+ r s _) 


y + r f z f 





(r 7 /2) • (r s+s+ /r fs+ r s+ - r s _ s _/r f - s _T s _) 





y + r f Zf 


r ff/ T ff 


Rm 


2rj, + rj(zfy — Wfx) 




r f 


r / ^ r f( z fy ~ w f x ) 


Fe±e±/T e+e - 


Rm 


3Rm 



The differences between the observed and the desired quantities are expected to be 0(1%). 

In principle, any D° branching fraction measured with DT yields in a coherent D°D° 
system is subject to such considerations. Analogous caveats pertain to B° and B Q S branching 
fractions measured with coherent B°B° and B®B° S pairs. However, for K°K° decays, such 
as those studied by KLOE, the situation is generally simpler. There, the desired branching 
fractions [5(J 5l[ are for Kg and Kl, rather than for K° and K°, so corrections for the 



lifetime asymmetry (y w 0.997) need not be applied. 



IV. EXPERIMENTAL SENSITIVITY 

The least-squares fit discussed in the previous section can be extended to extract the 
parameters y, x 2 , rf, Zf, and WfX (C-even only), in addition to Bj and Af. Efficiency- 
corrected ST and DT yields are identified with the functions given in Tables [TH IHIl and IIVI 
We estimate uncertainties on the fit parameters based on approximately 3 x 10 6 D°D° pairs, 
using efficiencies and background levels similar to those found at CLEO-c. In the rate 
ratios in Table |V] uncertainties that are correlated by final state, such as tracking efficiency 
uncertainties, cancel exactly. Therefore, the uncertainties on the mixing and amplitude ratio 
parameters stem primarily from statistics and from uncorrelated systematic uncertainties. 

The decay modes considered are listed in Table IVIl There exist, in principle, different 
rf, Zf, and Wf parameters for each mode / included in the fit. Therefore, for simplicity, we 
include only one such mode in the analysis: D — ► fC ± 7r T . In practice, adding more hadronic 
flavored modes does not noticeably improve the precision of y because the limiting statistical 
uncertainty comes from DT yields involving S±. The branching fraction determinations, 
however, would benefit from having additional modes in the fit. 

We omit ST yields for modes with a neutrino or a K\\, which typically escapes detection, 
because they are difficult to measure. In principle, one could reconstruct the remainder of 
the event inclusively to infer the presence of the missing particle from energy and momentum 
conservation j52|. The efficiency of this method depends on the hermeticity of the detector. 

For DT modes with one missing neutrino or K\, this method is more straightforward 
to implement because the D° and D° are both reconstucted exclusively. Therefore, we do 
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TABLE VI: Final states included in the fits, along with assumed branching fractions, signal ef- 
ficiencies, and expected single tag yields for N = 3 x 10 6 . ST yields in square brackets are not 
included in any of the fits. 



Final State 


Type 


B (%) 


e(%) 


ST Yield (10 3 ) 


K--K+ 


/ 


3.91 


66 


78 


K+TT- 


/ 


3.91 


66 


78 


K-e + v e 


e+ 


3.5 


62 


[65] 


K + e~u e 




3.5 


62 


[65] 


K + K~ 


s+ 


0.389 


59 


14 




s+ 


0.138 


73 


6.0 


l^vrV 


s+ 


0.89 


15 


8.0 


K° l tt 


s + 


1.15 


62 


[43] 


(K° S TT+7T-)CP+ 


s+ 


1.0 


38 


27 


(K° L 7r+TT-)cp + 


s+ 


1.0 


76 


[46] 


K° s <t> 




0.47 


7.7 


2.2 


K Q s uo 


s_ 


1.15 


14 


9.7 


K S 7T° 




1.15 


31 


21 


K° L TT Tr 


S- 


0.89 


30 


[16] 


(K° S TT+TT-) CP _ 


5_ 


1.0 


38 


23 


(K° L ir+ir-) C P- 


S- 


1.0 


76 


[46] 



include these yields in the fit. In the case of a missing K®, one must veto Kg decays, which 
would have the opposite CP eigenvalue. 

For DT modes with two undetected particles, one can constrain the event kinematically, 
up to a twofold ambiguity j53|. Background events tend to fail these constraints, so the 
signal can be isolated effectively. We assume this method is used to measure £ + £~, £ ± £ ± , 
KIP, and K° L K° L DT yields. 

The input ST yields are listed in Table IVIt and the input DT yields are derived from 
products of the ST branching fractions and efficiencies. In most cases, the background is 
negligible (cx^ at ~ \N). We also include a conservative 1% uncorrelated sytematic uncer- 
tainty on each yield measurement. For modes that only have contributions from Rm, which 
we assume to be zero, we use yield measurements of ± 1 ± 1. Yields for forbidden modes 
(S±S± for C-odd, S + S- for C-even) are not included. The fit accounts for the statistical 
correlations among ST and DT yields. 

We perform fits for both C eigenvalues using equal numbers of D°D° pairs. To improve 
the precision of these fits, we include external measurements of branching fractions 48j. The 



first two columns of Table IVHI show the expected uncertainties on the mixing and strong 
phase parameters from these fits. The dramatic difference between the y uncertainties for 
the two cases stems from the negative correlation between Bs + and £>s_ introduced by the 
presence of S + S_ yields for C-odd; these branching fractions are positively correlated for 
C-even. From the second line of Table EJ it can be seen that a negative correlation increases 
the uncertainty on y. The difference between the x 2 uncertainties for the two cases reflects 
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TABLE VII: Estimated uncertainties (statistical and systematic, respectively) for different C con- 
figurations, with branching fractions constrained to the world averages. We include C-even ST 
yields in the second column, but not the third. 



Parameter 


Value 


M°- = 3 x 10 6 


M c+ = 3 x 10 6 


M°- = 10 • N c+ = 3 x 10 6 


y 





±0.015 ± 0.008 


±0.007 ± 0.003 


±0.012 ± 0.005 


x 2 (1(T 3 ) 





±0.6 ±0.6 


±0.3 ±0.3 


±0.6 ±0.6 


COS 5ktt 


1 


±0.21 ±0.04 


±0.27 ±0.05 


±0.20 ± 0.04 


X Sin 5 Kn 







±0.022 ± 0.003 


±0.027 ± 0.005 


r 2 (10~ 3 ) 


3.74 


±1.0 ±0.0 


±1.7 ±0.1 


±1.0 ±0.0 



the factor of 3 accompanying Rm in F + (£ , ^ ± ). 

In reality, C-even D°D° pairs [from e + e~ — > D°D 0/ y(n7r )} are produced above D*D 
threshold and are more difficult to identify than C-odd D°D° pairs produced at the ^(3770). 
In particular, one must distinguish D°D ^( y mr ), which gives C-even, from D D (mr ) and 
_D _D 77, which give C-odd. While it is possible to do so for DT modes, where the entire 
event is reconstructed, a large uncertainty is incurred for ST modes. Therefore, we perform 
a third fit that combines yields from both C configurations, but with C-even ST yields 
omitted. Because of the smaller cross section and efficiencies for D°D°'~f (only half the soft 
photons can be identified), C-even DT yields above D*D threshold are an order of magnitude 
smaller than C-odd yields from an equal luminosity at the ?/>(3770). Results for this fit are 
also shown in Table IVHl 

One important source of systematic uncertainty not included above is the purity of the 
initial C = ±1 state. The sample composition can be determined from the ratios 

T(S + ,S' + )T(S-,S'_) _ T(S + ,S' + )T(S-,SL) _ 4T(S + , S + )T(S-, gj _ (Af c+ \ 2 

T{S+, S_)T{S' + , S'_) T(S+, S'_)T(S' + , SJ) T 2 (S + ,S_) \Af c ~ J ' 1 ' 

assuming CP is conserved. Thus, if we include the forbidden S+S+, SS-, and S+S- DT 
yields that were previously omitted, then we can construct every other ST or DT yield as a 
sum of C-odd and C-even contributions, with their relative sizes constrained by Equation [361 
In this way, the systematic uncertainty is absorbed into the statistical uncertainties, and the 
C content of the sample is self-calibrating. If the fit is performed on pure samples, with 
either T(S±,S±) or T(S + ,S-) measured to be consistent with zero, then the uncertainties 
given in Table IVIII suffer no degradation. On the other hand, if N c ~ = N c+ , then there 
is an unfortunate cancellation of terms containing Zf or y, and these parameters cannot be 
determined at all. 

To demonstrate the importance of semileptonic modes in this analysis, we consider a 
variation on the above fits. If no semileptonic yields are measured, then, for C-odd there is 
only one independent combination of 77, Zf, and?/: [r fZf+y{l+r 2 )]/ {l+r 2 +r fZfy) w rfZf+y. 
On the other hand, if only the i^K® DT modes are omitted, then y and Zf can be determined 
separately but with uncertainties approximately 50% larger than those in Table IVHl where 
the i^K® DT modes are included. 

The £ ± £ ± and £ + £~ DT modes improve the uncertainties on x 2 and Be, but not on y 
and Zf. However, if the self-calibrating fit described above were performed without these 
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modes, then rj would be strongly coupled to x 2 and Zf [through r(/, /)], resulting in inflated 
uncertainties. In order to stabilize the fit, it would be necessary to fix the value of x 2 . 

Experimental constraints on charm mixing are usually presented as a two-dimensional 
region either in the plane of x' Kn versus y' K7T or in the plane of x versus y. In Figures [I] 
and [21 we compare current constraints with those projected using the method described in 
this paper. In Figured], we show the results of the time-dependent analyses of D° — > K + tt~ 
from CLEO [31], BABAR Belle @, and FOCUS 0. The regions for CLEO, BABAR, 
and Belle allow for CP violation in the decay amplitude, in the mixing amplitude, and in 
the interference between these two processes, while the FOCUS result does not. Several 
experiments have also measured y directly by comparing the D° decay time for the K~n + 
final state to that for the CP eigenstates K + K~ and 7r + 7r~. The allowed region for y (labeled 



AT) shown in Figures Q] and [2] is the average of the results from E791 25|, CLEO [271 ]. 
BABAR 29], and Belle j30|. In depicting the y and Dalitz-plot results in Figure [U we 
assume 5k-k = 0; a non-zero value for Sktt would rotate the D° — > K~ K + /tt + 7t~ confidence 
region clockwise about the origin by an angle Sxn- The best limit on Rm from semileptonic 
searches for charm mixing (D° — > D° — > K + £~i>i), shown in Figures [T] and [2} is from 



the Belle experiment 
D° -» K + 7i~n° and D° -* K + n 
and CLEO 123] are not shown. 



241 ] . Figure [T] also displays the BABAR Rm limits from wrong-sign 



"7T + 7T" 



Semileptonic results from E791 [2l|, BABAR [22 



In Figure [2j we plot the projected 95% confidence level (C.L.) contour for the fit in the 



third column of Table IVIH along 
of D° -> K° s n + 7i- by CLEO ET 



with the results of a time-dependent Dalitz-plot analysis 
as well as the lifetime and semileptonic results discussed 
above. We note that our sensitivity to x depends strongly on the value of 8k-k, while our 
sensitivity to y does not, and lowering the value of | cos^,,-] would reduce the area of the 
contour. In general, our expected upper limits compare favorably with the current best 
limits on charm mixing. 



V. SUMMARY 

We have derived ST and DT rate expressions for correlated D°D° pairs in a definite 
C eigenstate. Interference between amplitudes for indistinguishable final states enhances 
some D° decays and suppresses others, depending on the D°-D° mixing parameters and on 
the magnitudes and phases of various amplitude ratios. By examining different types of 
final states, which have different interference characteristics, we can extract the mixing and 
amplitude ratio parameters, in addition to the branching fractions. In contrast to previous 
measurements of mixing parameters 0] , our method is both time- independent and sensitive 
to x and y at first order, so it is subject to different systematic uncertainties. Also, it 
is unique to threshold production where the D°D° initial state is known, unlike with D° 
mesons produced at fixed target experiments or through D* decays, as at the I?-factories 
and at LEP. 

When performing this analysis for K°K°, B°B°, and B®B° S decays, one should incorporate 
CP violation and non-trivial weak phases. In principle, given a source of coherent B Q S B° S 
pairs, such as those produced through the T(5S), B° s -B° s mixing could be probed in a fashion 
similar to that presented in this paper for D°-D° mixing. 
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FIG. 1: Allowed regions in the plane of y' Kn versus x' Kn |2J, |25j, ]21, M, 120, M, \M, \M, ISJ, m, \M, \m > 
assuming 5xn = 0. A non-zero value for 5k-k would rotate the Ar confidence region clockwise about 
the origin by an angle 8kk- 
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